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By molecular dynamics simulation, we have investigated classical Heisenberg spins, which 
are arrayed on a finite simple cubic lattice and interact with each other only by the dipole- 
dipole interaction, and have found its peculiar from- Edge-to -interior freezing process. As the 
temperature is decreased, spins on each edge predominantly start to freeze in a ferromagnetic 
alignment parallel to the edge around the corresponding bulk transition temperature, then 
from each edges grow domains with short-range orders similar to the corresponding bulk 
orders, and the system ends up with a unique multi-domain ground state at the lowest 
temperature. We interpret this freezing characteristics is attributed to the anisotropic and 
long-range nature of the dipole-dipole interaction combined with a finite-size effect. 

KEYWORDS: Magnetic dipolar system, nanomagnetism, magnetic domain, spin dynamics sim- 
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In the last decade systems consisting of arrayed ferromagnetic nanoparticles have attracted 
much attention as a possible element with the high storage density. 1 Magnetic properties of 
such systems have been analyzed based on various theoretical model, in which the magnetic 
moment of each nanoparticle is represented by a classical Heisenberg spin with a proper 
magnitude. The spins are interacting with each other by the dipole-dipole interaction, and 
each spin suffers from the magnetic anisotropy energy which represents the shape and/or bulk 
lattice anisotropics of the original nanoparticle magnetic moments. There appeared studies on 
the susceptibility, 2 magnetic hysteresis, 3 ' 4 and energy relaxation 5 in such models of a finite 
size. However, there have been little work on systems only with the dipole-dipole interaction 
which we call here simply dipolar systems. Considering that to clarify magnetic properties 
of dipolar systems of a finite size is of importance as the first step in researches of arrayed 
magnetic nanoparticles as well as of interest from a theoretical viewpoint in statistical physics, 
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Fig. 1. The spin configurations in the ground state, (a) Antiferromagnetically aligned ferromagnetic 
chains (af-FMC) order, and (b) one of the LT orders which have the continuous 0(3) symmetry. 

we address this problem in the present work. 

The ground-state magnetic structure of the bulk dipolar system was studied more than 
a half century ago. 6 ' 7 In particular, Luttinger and Tisza (LT) 7 derived analytically that of 
cubic lattices, restricting themselves to magnetic superlattice structures of 2x2x2 unit. The 
LT ground state on the simple cubic (SC) lattice is given by any properly-normalized linear 
combination of the three basic states (Fig. 1(a)), one of which is shown in Fig. 1(b): spins on 
z-axis align ferromagnetically along the axis and these ferromagnetic chains align antiferro- 
magnetically in the x and y directions (the rest two are defined by changing the directions 
x, z). We call them the antiferromagnetically aligned ferromagnetic chains (af-FMC) orders. 
Surprisingly, a very limited number of works on the model have been done since then. 8 ' 9 Re- 
cently we have checked the LT ground states by the Monte Cairo (MC) analysis. 10 We have 
found the lower energy states with the longer magnetic superlattice structures than those of 
the LT ones for the body centered cubic lattice, while for the SC lattice no lower energy state 
than the LT one has been found. 

In the present work we focus on a dipolar cube cut out from the SC dipolar lattice, 
with its edges being parallel to the lattice axes, and look for finite-size effects on its ground- 
state structures and its freezing characteristics by the stochastic molecular dynamics (MD) 
simulation. We have found very peculiar aspects which one cannot think of for finite-size 
systems with any short-range interaction. The ground state of the dipolar cube consists of 12 
domains, each of which is very close to one of three af-FMC orders mentioned above. A very 
peculiar aspect we have observed in the cooling process is what we call from-edge-to-interior 
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freezing characteristics: spins on the edges first freeze independently from spin freezing on 
other edges around the transition temperature of the corresponding bulk SC dipolar system, 
the af-FMC domains grow from the edges to the interior at lower temperatures, and then 
finally a unique combination of the domains is reached at the lowest temperature. 

The model we study consists of classical Heisenberg spins, S"s, which are arrayed on a 
finite cube. They correspond to magnetic moments of single ferromagnetic domain particles 
normalized by their magnetization M s . The dipole-dipole interaction energy between the spins 
is written as 

Here denotes dimensionless length (normalized by the lattice constant, a) between sites 
i and jf, and the unit vector along the direction from site i to site j. To the coupling 
constant, J, the value of M 2 is contained, and J /a 3 is put the energy (and temperature T 
with /cb = 1) unit of the present model. For example, in Ni nanoparticles array with 30 nm 
particle diameter and a ~ 100 nm, 3 ' J /a 3 is estimated about 5 K. If Co is replaced for Ni, 
J /a 3 becomes of order of several tens K. It is also worth quoting a typical value of magnetic 
molecular crystal, J /a 3 ~ 0.7 K for S z — 10 with a ~ 0.7 nm. 11 ' 12 

The freezing characteristics of the model is analyzed by the Landau-Lifshiz-Gilbert (LLG) 



eq., 13 



with the effective magnetic field, i? e ff,z5 given by 

^eff,z + fi; (3) 

oSi 

where H is the dipole-dipole interaction of Eq. (1). To include the heat-bath effect of temper- 
ature T we introduce random force its root mean square value being proportional to Ta. 
a and 7 are so-called Gilbert damping constant and the gyromagnetic constant, respectively. 
In the present study we solve the above set of equations based on the Euler method with 
conserving the norm of spins (M s = 1) and by setting aa 3 /M s ~fJ as the time unit. We further 
put both constants a and 7 to 1.0, and the time step of integrating Eq. (2), At, to 0.01. With 
this choice of the parameter values, the Larmor precession of a spin due to the internal field of 
averaged magnitude damps in a time comparable with one period of the precession at T = 0. 

In the present work we mainly report the results obtained by simulations on a cube whose 
linear size, L, is set to 16. It is an optimal size to look for peculiar phenomena in a magnetic 
dipolar cube: the size should be large enough to enable us to observe typical behavior of spins 
both on the boundary (corners, edges, and surfaces) and inside the cube, while simulation 
on a too large system is beyond computational power available for us. We note that, in the 
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present simulation the sum in Eq. (1) are faithfully taken into account except for the analysis 
to obtain the result shown in Fig. 3(b) below. 

We perform typically 8 cooling runs with a fixed cooling rate, in which we use different 
sets of random numbers generating their initial configurations as well as fluctuation forces. 
The average over them is carried out to obtain nearly-equilibrium, or frozen values of the 
observables reported below unless otherwise indicated. In each cooling run, temperature T is 
initially set to 1 (= J /a 3 ) and is decreased by a step of AT = 0.02. At each temperature 
Eq. (2) is integrated by a period of r — 5 x 10 4 At, which corresponds to about 500 periods 
of the Larmor precession mensioned above. The average value of quantities of our interest at 
each temperature is defined by that over the whole period of r. For example, the freezing of 
individual spins, is evaluated as 

(4) 

Let us first present the results of global thermodynamic quantities. In Fig. 2, we demon- 
strate the specific heat, C, and the overall freezing parameter, <S, as a function of temperature. 
Here S is defined as S = TV -1 X^<%> including the average over the cooling runs, where N is 
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the total number of spins in the cube. It represents on average how much each spin is frozen in 
a time scale of r. We see from Fig. 2 that the specific heat exhibits a peak at T* ~ 0.6, which 
is closed to the transition temperature, T c = 0.56, evaluated previously in the corresponding 
bulk system. 9 Correspondingly, S starts to rapidly increase around T*. Both C and S clearly 
indicate occurrence of cooperative freezing around T*. 

It is natural for us to expect that the frozen order we found in the cube somehow reflects the 
magnetic structure of the bulk ground-state derived analytically by LT. The order parameter, 
which represents the LT ground-state order (see Fig. l(a,b)), is defined by a vector, Mlt, 
whose a-component is written as 



where a, (3 G and nf denotes the ^-coordinate of site i. The absolute magnitudes 



rapidly increase from around T* with decreasing temperature and reach almost a unique 
value of about 0.8 which is significantly less than the value of the bulk order, 1.0, at the 
lowest temperature. One plausible interpretation of the result is that spins on the cube reach 
a stable state which consists of the domains with different af-FMC orders. Figure 3(b), on the 
other hand, are .A/fur's obtained by the identical simulation as above but with one difference, 
namely, the dipole-dipole interaction is cut off within the range of L/2 in this simulation. By 
the introduction of the cut-off, the freezing behavior of Aiur fluctuates drastically among 8 
cooling runs, though C and S do not exhibit such a qualitative difference (not shown). 

To check the above interpretation on the result shown in Fig. 3(a), we examine freezing 
patterns {Si} defined by Eq. (4) of whole the cube. Figures 4(a-d) are typical result observed 
in one of the cooling processes. At T = 0.68 (> T*), spins on some edges start to freeze. As the 
temperature further decreases (T = 0.6 < T*), the freezing extends to the whole boundary 
in the sense that boundary spins are relatively less fluctuating than inner spins. This feature 
is kept at further lower temperatures. At T = 0.52, however, spins in a domain indicated by 
the circle in Fig. 4(c) look as if they are partially melting. Finally at very low temperatures 
such as T = 0.36, also inner spins almost freeze in the time scale of our observation. We can 
say that, in a cooling process of the present dipolar cube, spins on edges are the most stable, 
and those on boundary surfaces are stabler than those inside the cube. 

A typical spin configuration at a very low temperature reached after the cooling process 
is shown in Figs. 5(a,b). A common property shared with such spin configurations is that 
spins on a boundary surface (including an edge) tend to align in the direction parallel to 
the surface (the edge). This tendency can be naively attributed to the anisotropic nature of 
the dipole-dipole interaction of Eq. (1). Namely, on a surface where there are two (one) axes 
parallel (perpendicular) to the surface, a ferromagnetic chain (FMC) aligned in parallel to 




(5) 



of TWlt, denoted as Mlt, for the 8 cooling runs are shown in Fig. 3(a). All TWlt's start to 
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the surface is more favorable. Naturally, on an edge an FMC order parallel to the edge is 
most favorable. Consequently, we can also expect that the magnetic structure at the lowest 
temperatures consists of 12 domains with short-range order which is close to one of the three 
af-FMC orders whose FMC order develops along the edge of each domain. This is in fact seen 
in Figs. 5(a,b). Also, if the simulation is started from one of the af-FMC configurations such 
as show in Fig. 1(a), it ends up with a similar configuration to the one shown in Figs. 5 (a, b): 
the spins on a boundary initially in the direction perpendicular to the boundary flop into the 
parallel directions to the boundary. 

By further inspection of the spin configurations, we find that the one with the lowest 
energy satisfies the following two conditions; the FMC orders on three edges which cross at 
each corner take 2-in(out)-l-out(in) configuration and those on opposite sides of a surface is 
aligned antiparallelly. These two conditions make the combination of domains in the ground 
state to be unique one which is shown in Fig 5(c). Here structures which can be connected by 
7r/2-rotation around 4- fold axis of the cubic sample are regarded as identical. In this context, 
we want to emphasize that the proper self-rearrangement of the domains is due to the very 
long-range nature of the dipole-dipole interaction. In fact, when the cut-off of the interaction is 
introduced, each cooling run ends up with states with different values of Aiur as indicated in 
Fig. 3(b), i.e., the adjustment of the af-FMC domains to reach their ground-state combination 
hardly takes place when the dipole-dipole interaction is cut off. 

One more important freezing aspect we have observed is as follows. In the ground-state 
configuration, individual local energies of the corner spins are higher than those of other 
spins as expected. However, those of the edge spins are not necessarily lower than those of 
the interior spins. This implies that the relatively stronger freezing tendency discussed above 
can not be simply attributed to energetics of individual spins. Instead, we have to consider 
cooperative nature of the freezing. As noted above, the bulk LT order has the continuous 0(3) 
symmetry, and so different LT short-range orders are expected to grow in the bulk system, 
or inside the cube as well, at higher temperatures than T c or T* (^ T c ). On the boundary of 
the cube, however, a symmetry breaking due to the presence of boundary partially solves the 
degeneracy. Especially on an edge, an FMC order along the edge direction, which has only 
the Z2 symmetry as a whole, becomes preferable. This symmetry lowering of the short-range 
order is considered to be an origin of the from- edge-to-interior freezing characteristics we have 
found at and below T*. Thus the boundary geometry, combined with the anisotropic nature 
of the dipole-dipole interaction, also plays an important role on the freezing nature of the 
dipolar cube. 

So far the dipolar cube with L = 16 has been discussed. For an odd-L dipolar cube, one 
of the two conditions for the ground-state configuration for an even-L cube mentioned above 
is certainly violated, since FMC's on opposite sides of an odd-L square in its ground state are 
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checked to be parallel to each other. 15 We also note the stable configurations in cubes with 
small L's. In the L = 16 cube, spins around a vertex of the cube take an optimized direction, 
which is regarded as a part of the so-called vortex structures. 8 In small cubes with L = 4, 6, 8, 
this corner structure extends to a size comparable with L/2, and consequently their ground 
state consists of a single-domain vortex state as shown in Fig. 3(d,e) which differs distinctly 
from the configuration shown in Figs. 5(a,b). We can expect variety of stable configurations 
of SC dipolar cubes arising from the relative angles between edges of the cube and the lattice 
axes. Further variety of magnetic properties of finite dipolar systems is expected when we 
think of the shape effect, different lattice structures, and so on. The present work has just 
opened a door in this interesting field. 

To conclude let us summarize a scenario of the from- edge-to-interior freezing process we 
have found in the dipolar cube. 1) At higher temperatures than T* (~ T c ), the LT short-range 
orders are expected to grow. 2) Around T*, where some of the LT short-range orders become 
in touch with one of the cube edges, the degeneracy in the LT orders is lifted, and one of 
the af-FMC orders fitted to the edge direction becomes stabilized. 3) The domains with the 
af-FMC order extends from edge to interior at lower temperatures. 4) In case the combination 
of domain structures does not match with the one of the ground-state, upturns of certain 
domains occur, and 5) the spin configuration reaches to the unique ground-state at lowest 
temperatures. 

The present work is supported by NAREGI Nanoscience Project from the Ministry of 
Education, Culture, Sports, Science, and Technology. The numerical simulations have been 
partially performed also by using the facilities at the Supercomputer Center, Institute for 
Solid State Physics, the University of Tokyo. 
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Fig. 3. The LT order parameter A^lt of the dipolar cube observed in 8 freezing runs; a) without and 
b) with cut-off for the dipole-dipole interaction. 
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Fig. 4. Patterns of {5*}'s defined by eq. 4 for T=0.625, 0.5, 0.375 and 0.125. 
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Fig. 5. The spin configurations for the surfaces (a) and interior (b) of the dipolar cube at lowest 
temperatures. The shaded regions denote domain walls. Figure (c) represents the domain order in 
the ground state of the dipolar cube. Figures (d,e) are the spin configuration in the ground state 
of the L=4 dipolar cube. 
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